Large Scale Pressure Fluctuations and Sunyaev-Zel'dovich Effect 



O 

o 
o 

(N 



in 



Asantha Cooray 

Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637. 
E-mail: asante Qhyde. uchicago. edu 
(Submitted to PRD - May 2000) 

The Sunyaev-Zel'dovich (SZ) effect associated with pressure fluctuations of the large scale struc- 
ture gas distribution will be probed with current and upcoming wide-field small angular scale cosmic 
microwave background experiments. We study the generation of pressure fluctuations by baryons 
which are present in virialized dark matter halos, with overdensities > 200 and by baryons present in 
overdensities < 10. For collapsed halos, assuming the gas distribution is in hydrostatic equilibrium 
with matter density distribution, we predict the pressure power spectrum and bispectrum associated 
with the large scale structure gas distribution by extending the dark matter halo approach which 
describes the density field in terms of correlations between and within halos. The projected pressure 
power spectrum allows a determination of the resulting SZ power spectrum due to virialized struc- 
tures. The unshocked photoionized baryons present in smaller overdensities trace the Jeans-scale 
smoothed dark matter distribution. They provide a lower limit to the SZ effect due to large scale 
structure in the absence of massive collapsed halos. We extend our calculations to discuss higher 
order statistics, such as bispectrum and skewness in SZ data. The SZ-weak lensing cross-correlation 
is suggested as a probe of correlations between dark matter and baryon density fields, while the 
probability distribution functions of peak statistics of SZ halos in wide field CMB data can be used 
as a probe of cosmology and non-Gaussian evolution of large scale pressure fluctuations. 



> 
oo 

(N 

in ■ 

o : 
o . 
o ■ 

^ ' 

6 ■ 
a: 



I. INTRODUCTION 

In recent years, increasing attention has been given 
to the physical properties of the intergalactic warm and 
hot plasma gas distribution associated with large scale 
structure and the possibility of its detection (e.g., 0). 
It is now widely believed that at least ~ 50% of the 
present day baryons, when compared to the total baryon 
density through big bang nucleosynthesis, are present in 
this warm gas distribution and have remained undetected 
given its nature (e.g., Q). Currently proposed methods 
for the detection of this gas with include observations 
of the thermal diffuse X-ray emission (e.g., ||), associ- 
ated X-ray and UV absorption and emission lines (e.g., 
|i|) and resulting Sunyaev-Zel'dovich (SZ; |J) effect (e.g., 

II). 

The SZ effect arises from the inverse-Compton scat- 
tering of CMB photons by hot electrons along the line of 
sight. This effect has now been directly imaged towards 
massive galaxy clusters (e.g., where temperature of 

the scattering medium can reach as high as 10 keV pro- 
ducing temperature changes in the CMB of order 1 mK 
at Rayleigh- Jeans wavelengths. Previous analytical pre- 
dictions of the resulting SZ effect due to large scale struc- 
ture have been based on either through a Press-Schechter 
(PS; M) description of the contributing galaxy clusters 
(e.g., ]lO| , pd| ) or using a biased description of the pressure 
power spectrum with respect to the dark matter den- 
sity field (e.g., ||). Numerical simulations (e.g., |l2 14|) 
are beginning to improve some of these analytical predic- 
tions, but are still limited to handful of simulations with 
limited dynamical range and resolution. Therefore, it is 
important that one consider improving analytical models 
of the large scale structure SZ effect, and provide predic- 



tions which can be easily tested through simulations. 

Our present study on the large scale baryon pressure 
and the resulting SZ effect is timely for several reasons, 
including the fact that improving numerical simulations 
have recently begun to make detailed predictions for the 
pressure power spectrum and SZ effect such that those 
predictions can be extended and improved with analyti- 
cal models ]l2|-|l4|]. Also, several studied have considered 
the possibility that large scale baryon distribution can 
be probed with upcoming CMB missions using SZ effect 
(e.g., H). Our calculations can be used to further refine 
these predictions and to investigate the possibility how 
such analytical model as the one presented here can be 
tested with observations. 

As part of this study, we extend previous studies by 
considering the full power spectrum and bispectrum, the 
Fourier space analog of the three-point function, of pres- 
sure fluctuations. The pressure power spectrum and 
bispectrum contains all necessary information on the 
large scale distribution of temperature weighted baryons, 
whereas, the SZ power spectrum is only a projected mea- 
surement of the pressure power spectrum. This can be 
compared to weak gravitational lensing, where lensing 
is a direct probe of the projected dark matter density 
distribution. The bispectrum of pressure fluctuations, 
and SZ bispectrum, contains all the information present 
at the three-point level, whereas conventional statistics, 
such as skewness, do not. An useful advantage of using 
the 3d statistics, such as the pressure power spectrum, 
is that they can directly compared to numerical simula- 
tions, while only 2d statistics, such as the projected pres- 
sure power spectrum along the line of sight, basically the 
SZ power spectrum, can be observed. Our approach here 
is to consider both such that our calculations can eventu- 
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ally be compared to both simulations and observations. 

The calculation of pressure power spectrum and bis- 
pectrum requires detailed knowledge on the baryon dis- 
tribution, which can eventually be obtained numerically 
through hydrodynamical simulations. Here, we provide 
an analytical technique to obtain the pressure power 
spectrum and bispectrum by describing the baryon dis- 
tribution in the universe as (1) present in virialized halos 
with overdensities > 200 with respect to background den- 
sities (2) unshocked diffuse baryons in overdensities < 10 
that trace a Jeans-smoothed dark matter density field 
(3) the intermediate overdensity region, which is likely 
to be currently undergoing in shock heating and falling 
on to structures such as filaments. In the present paper 
we discuss the first two regimes, while a useful approach 
to include the latter, through simulations, is discussed. 

Our description of baryons present in virialized halos 
follow recent studies on the dark matter density field 
through halo contributions [fl5| |l7| following [|l8| and ap- 
plied to lensing statistics in |19[ and |20[ . For the de- 
scription of baryons, the critical ingredients are: the PS 
formalism j|] for the mass function; the NFW profile of 
PH , and the halo bias model of p2fl . The baryons are 
assumed to be in hydrostatic equilibrium with respect to 
dark matter distribution, which is a valid assumption, 
at least for the high mass halos that have been observed 
with X-ray instruments, given the existence of regularity 
relations between cluster baryon and dark matter phys- 
ical properties (e.g., |2j|). We take two descriptions of 
the temperature structure: (1) virial temperature and (2) 
virial temperature plus an additional source of minimum 
energy. The latter consideration allows the possibility for 
a secondary source of energy for baryons, such as due to 
preheating through stellar formation and feedback pro- 
cesses. Numerical simulations (e.g., pp^]), as well ob- 
servations (e.g. , |25|,|26| ) , suggest the existence of such an 
energy source. The low photoionized overdensity baryons 
are described following the analytical description of [^7) . 
The fraction of baryons present in such low overdensities 
are assumed to follow what has been measured in nu- 
merical simulations of JlJ . We suggest that such baryons 
provide a lower limit to the SZ effect in the absence of 
any contribution from baryons present in virialized halos. 

Throughout this paper, we will take ACDM as our fidu- 
cial cosmology with parameters f2 c = 0.30 for the CDM 
density, = 0.05 for the baryon density, £!a = 0.65 for 
the cosmological constant, h = 0.65 for the dimensionless 
Hubble constant and a scale invariant spectrum of pri- 
mordial fluctuations, normalized to galaxy cluster abun- 
dances (as = 0.9 see ^8|) and consistent with COBE 
p9fl . For the linear power spectrum, we take the fitting 
formula for the transfer function given in |3C|| . 

The paper is organized as following: In §[n], we re- 
view the dark matter halo approach to modeling the den- 
sity field and extend it to model properties associated 
with large scale baryon distribution, mainly the pressure 
fluctuations that contributes to the observable SZ effect. 
We suggest recent papers by Seljak p5|, Ma & Fry p6[, 



Cooray & Hu PQ| , and Scoccimarro et al |17]] for details 
on the dark matter halo approach and applications to 
other observable statistics such as galaxy properties and 
weak gravitational lensing. As necessary, we use tech- 
niques developped in these papers for our current calcu- 
lation. In §[Hl| we apply the formalism to the convergence 
power spectrum, skewness, and bispectrum. We conclude 
in §IV with a summary of our main results. 



II. DENSITY AND PRESSURE POWER 
SPECTRA 

A. General Definitions 

In order to calculate the contribution to temperature 
anisotropies through SZ effect associated with large scale 
structure, we divide the LSS with overdensities > 200 
as collapsed and virialized halos with a gas distribution 
following hydrostatic equilibrium and with virial temper- 
atures. The pressure power spectrum can be calculated 
using an extension to the dark matter halo approach by 
assuming a physical relation between baryons and dark 
matter. 

The baryons with overdensities < 10 track the dark 
matter distribution and their power spectrum has been 
studied by [ p7fl . These baryons have temperatures sim- 
ilar to photoionization energies of Hydrogen and He- 
lium. Current numerical simulations, e.g., Q, suggest 
that most of the baryons are in such low overdensities at 
z > 1, while at present day, are in virialized halos. The 
calculation of SZ effect due to such baryons follow ||, 
except that we include the redshift dependence of mass 
fraction within such low density halos following the nu- 
merical results of |lj , and modify the mean temperature 
of such baryons to be consistent with photoionization en- 
ergies. 

First we discuss the pressure and related power spectra 
due to collapsed halos. 

The dark matter profile of collapsed halos are taken to 
be the NFW pl| with a density distribution 



Ps(r) 



{r/r s )(l+r/r s y 



(1) 



The density profile can be integrated and related to the 
total dark matter mass of the halo within r v 



M s = 4np s r 3 s 



log(l + c) - 



1 



(2) 



where the concentration, c, is r v /r s . Choosing r v as the 
virial radius of the halo, spherical collapse tells us that 
M = 47rr^A(z)pfc/3, where A(z) is the overdensity of 
collapse (see e.g. fl3lf ) and pb is the background matter 
density today. We use comoving coordinates throughout. 
By equating these two expressions, one can eliminate p s 
and describe the halo by its mass M and concentration 
c. 
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FIG. 1. Pressure power spectra with single halo (PP; dotted line), double halo (hh; dashed line) and total (t; solid line) 
contributions, (a) Here, we use virial temperature to de scribe electrons and show pressure bias, b 7r (k) (long dashed line), filtered 
non-linear dark matter density power spectrum (§ II G ) and the measured pressure power spectrum in simulations (thick solid 
line), (b) The variations in pressure power with maximum mass. The dot-dashed lines show the total power with a minimum 
electron temperature of 0.75 keV, as an attempt to reproduce power spectra under a possible preheating scenario. 



Following |li| , we take the concentration of dark mat- 
ter halos to be 



c(M,z) = a(z) 



M 



-6(2) 



(3) 



where a(z) = 10.3(1 + z)-°- 3 and b{z) = 0.24(1 + z)-°- 3 . 
Here M*(z) is the non-linear mass scale at which the 
peak-height threshold, v(M,z) — 1. The above con- 
centration is chosen so that dark matter halos provide 
a reasonable match to the the non-linear density power 
spectrum as predicted by the it extends the treat- 
ment of |ll| to the redshifts of interest for SZ effect. We 
caution the reader that eqn. (||) is only a good fit for the 
ACDM model assumed. 

The gas density profile, p g (r), is calculated assuming 
the hydrostatic equilibrium between the gas distribution 
and the dark matter density field with in a halo. This 
is a valid assumption given that current observations of 
halos, mainly galaxy clusters, suggest the existence of 
regularity relations, such as size-temperature (e.g., p3[), 
between physical properties of dark matter and baryon 
distributions. 

The hydrostatic equilibrium implies, 



kT e d\ogp g GM s (r) 



finip 



dr 



(4) 



where now the M$(r) is the mass only out to a radius 
of r. Note that we have assumed here an isothermal 
temperature for the gas distribution. Solving for the the 
equations above, we can analytically calculate the baryon 
density profile p g (r) 



Pg( r ) = Pgoe 



-6 



br s /r 



where 6 is a constant, for a given mass, 



(5) 



(6) 



with the Boltzmann constant, ks- 

In general, the halos are described with virial temper- 
atures 



k B T e = 



~fGpm p Ms(r v ) 
3r v 



with 7 = 3/2 and /i = 0.59. In addition, we also con- 
sider the possibility for the existence of a constant non- 
gravitational energy in small mass halos, consistent with 
observations of galaxy groups, due to what is commonly 
known as "preheating" . The possibility for such a min- 
imum energy for baryons today comes from heating be- 
fore virialization due to energy injection and feedback 
processes such as processes associated with stellar for- 
mation. The total gas mass present in a dark matter 
halo within r v is 



M g (r v ) = 47rp g0 e 



dxx 2 (l + x) b/x . 



(8) 



The physical properties of the profile defined in Eq. (||) 
for baryons within dark matter halos, and a comparison 
to commonly used profiles such as isothermal and so- 
called beta-profiles, can be found in ||| and pi] . Com- 
pared to conventional profiles, this profile has the advan- 
tage that it is directly related to the dark matter profile 
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parameters, such as central density p s and concentration 
via scale radius r s , thus, any changes to the dark mat- 
ter distribution produces resulting changes in the baryon 
distribution. Also, one can study the effect of tempera- 
ture variations on the gas distribution as the parameter 
b defined in Eq. (|^) depends on it. A proper normal- 
ization for the dark matter halo distribution containing 
baryons comes through the c(M, z) relation in Eq. (|3j) 
such that the non-linear dark matter power spectrum is 
produced in numerical simulations by same halos. Note 
that our hydrostatic equilibrium ignores the self-gravity 
contribution from baryons to the total potential as we 
only include the dark matter contribution to total mass. 
Since the baryon mass is expected to be < 10% of the 
total mass, we can safely ignore, as a first approximation, 
the contribution to total mass from baryons themselves. 

Roughly speaking, the perturbative aspect of the clus- 
tering of the dark matter and baryons is described by the 
correlations between halos, whereas the nonlinear aspect 
is described by the correlations within halos, i.e. the halo 
profiles. We will consider the Fourier analogies of the two 
and three point correlations of the dark matter density, 
<5, baryon pressure, II, and galaxy distribution, g, defined 
in the usual way 



(^(k)^(k')) = (2^(k-k')^(fc), 



(9) 



(&(ki)&(k2)fc(k3)) = (2 7 r) 3 (5(k 1 + k 2 + k 3 )B*(fci, fe, k 3 ) . 

(10) 

with i representing S, II or g. We will also consider cross- 
correlations between the two, such as the dark matter 
density-pressure power spectrum Pg U (k), which is what 
one probes by correlating, say, the SZ effect and weak 
gravitational lensing observations. Here and throughout, 
we occasionally suppress the redshift dependence where 
no confusion will arise. 

As presented in ^0|, these spectra are related to the 
linear density power spectrum P(k) through the bias pa- 
rameters and the normalized 3d Fourier transform of the 
density profile Pi(r, M) 



dr 4nr 2 pi(r, M) 



sin(fcr) 
kr 



(11) 



where i represents either the density, <5, or the gas, g, 
profile and associated masses respectively given in equa- 
tions and ^| With an increase in temperature rel- 
ative to the virial temperature of the halo, especially 
for halos with masses < 10 13 Mq, the gas profile is 
such that it does not fall rapidly at the virial radius, 
leading to an arbitrary cut off when doing the Fourier 
transformation. We included an additional filter to the 
gas density profile such that the gas density profile de- 
creases smoothly but promptly to zero at the virial ra- 
dius: p' (r) = p(r)[erfc(r — r. u /\/2Ar) — 1] with Ar < r s . 
Detailed aspects of the pressure and SZ statistics due to 



medium to small mass halos (< 10 13 Mq) are sensitive 
to sharpness of this transition, but these issues do not 
change our primary results. Here, we concentrate mostly 
on the statistics due to massive and rare halos. Another 
possibility not considered here is to include the role of 
baryons at the outskirts of halos. The physical properties 
of such baryons can be semi-analytically calculated fol- 
lowing the assumption that virial radius provides a shock 
boundary for the equilibrium of baryons within and out- 
side virialized regimes. The baryons outside halos are 
likely to preheated and trace the Jeans-smoothed version 
of the dark matter density field. The proper inclusion of 
such baryons require the aid of numerical simulations or 
semi- analytical models. These baryons are likely to in- 
clude the ones present in overdensities between 200 and 
10, which we have neglected in the present calculation. 

Following po|| , it is convenient to define a general in- 
tegral over the halo mass function dn/dM , 



M n„ 



x T e (M, z) v y il (ki, M) . . . ytAk^M) , (12) 



where &o = 1. (2^| gives the following analytic predictions 
for the bias parameters which agree well with simulations: 



h(M;z) = 1 + 



v 2 {M;z) - 1 

So 



and 



v 2 (M-z) - 3 



(13) 



(14) 



Here, T e (M, z) is the electron temperature of the baryon 
distribution of given halo when pressure power spectrum 
is considered, v(M, z) — 5 c /a(M, z), where a(M,z) is 
the rms fluctuation within a top-hat filter at the virial 
radius corresponding to mass M, and S c is the threshold 
overdensity of spherical collapse (see |3l]] for useful fitting 
functions) . 

The terms related to the galaxy power spectrum in- 
cludes the average number of galaxies per dark matter 
halo, (N g ), the mean number density of galaxies in the 
universe n g . These parameters are discussed in § HE 
involving the galaxy-pressure power spectrum. 

We use the Press-Schechter (PS; mass function to 
describe dn/dM. We take the minimum mass to be 10 3 
Mq while the maximum mass is varied to study the effect 
of massive halos on related statistics. In general, masses 
above 10 16 Mq do not contribute to low order statistics 
due to the exponential decrease in the number density of 
such massive halos. 

To summarize, in comparison to previous work on the 
SZ effect from virialized halos, our model has following 
advantages: 

(1) A physically motivated profile for the distribution of 
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FIG. 2. Pressure-density cross-correlation power spectra, (a) Comparison of power spectra, rn(fc) is the correlation between 
dark matter and baryon distributions, (b) The variations in pressure power by changing the maximum mass considered. The 
figure follows Fig. ty. 



baryons in virialized dark matter halos, instead of an as- 
sumed profile such as the isothermal model or so-called 
/3-profile. 

(2) A mass function for the virialized structures with the 
dark matter distribution of such halos reproducing the 
numerically simulated dark matter power spectrum and 
higher order correlations. 

(3) Easily extendable variations to the baryon physics so 
as to account for issues such as pre-heating. 

(4) Direct calculation of 3d properties of the large scale 
baryon distribution, such as the pressure power spectrum 
and bispectrum, which can be compared easily in numer- 
ical simulations. 

We now discuss the calculation of properties related to 
the dark matter and baryons. In the case of baryons, we 
discuss pressure as this is the property that leads to the 
SZ effect allowing a useful probe of them. 



B. Density Power Spectrum 

Following [jl5j, we can decompose the density power 
spectrum, as a function of rcdshift, into contributions 
from single halos (shot noise or "Poisson" contributions), 



and correlations between two halos 



7-0,0,0/ 



2,55 ■ 



phh 
^55 



(A;) 



Ilf'°(k) P(k), 



such that 



pt _ p pp , phh 

~XX — r XX T r 88 ■ 



As k 



88 — 1 88 
P(k). 



(15) 



(16) 



(17) 



C. Pressure Power Spectrum 

As above, we can decompose the pressure power spec- 
trum, as relevant for the SZ effect, into contributions 
from single halos 



^nn(k) — I'llggikj k) , 
and correlations between halos 



p/i/i 
r nn 



(k) 



J#°(fc) P(k) 



such that 



Pnu(ky = P^(k) + P^(k) 



hh i 



(18) 



(19) 



(20) 



For the pressure power spectrum, since r\ in Eq. ([12]) is 
non-zero, there is additional mass weighing arising from 
the fact that T e oc M 2 / 3 resulting in an additional mass 
dependence. The dependence is such that most of the 
contributions to the pressure, and thus, to the SZ, power 
spectrum comes from most massive and rarest halos. 
This dependence has already been observed in numeri- 
cal simulations by [[uj. 

As we discuss later, this dependance on high mass ha- 
los to produce most of the pressure fluctuations also leads 
to several interesting results with regards to the detec- 
tion and observability of SZ effect, among which are 

(1) Most of the contribution to large scale SZ effect re- 
sults from massive clusters of galaxies, while smaller mass 
halos and structures at low electron temperatures such as 
filaments do not contribute significantly 

(2) Since massive halos dominate the SZ effect, and the 
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distribution of these halos are Poisson and highly non- 
Gaussian, most of the contribution to two point and 
higher-order statistics of SZ effect will be dominated 
by Poisson term and there will be a significant non- 
Gaussianity associated with large scale SZ effect 
In fact, as we find later, the large scale correlations only 
contribute at a level of 10% to the SZ power spectrum 
suggesting that such correlations can be mostly disre- 
garded. The non-Gaussianity associated with SZ effect 
may become a useful tool to separate out its contribu- 
tion from other sources of foregrounds in CMB anisotropy 
data, though this task can be efficiently carried out us- 
ing frequency information (see, J(|). We will return to all 
these issues in later sections. 



E. Galaxy-Pressure Power Spectrum 

The cross-correlation between the galaxy distribu- 
tion and gas field, as appropriate for galaxy-SZ cross- 
correlation can be decomposed as a single halo 



p pp 



(k) 



h',gS (k, fc) ' 



and 



jjhh 



(k) 



i$\k) p(k), 



such that 



Pns(ky=p£f(k)+P^(k) 



(28) 



(29) 



(30) 



D. Density-Pressure Power Spectrum 

The cross-correlation between the density and gas field, 
as appropriate for lensing-SZ cross-correlation can be de- 
composed as a single halo 

tjPP/,\ J-0,1,0/ 



(k)=I^(k,k) 



and 

Pns(k) = 
such that 

Pas{k)* 



r l,0,0/ 
' 1,5 



(fc) P(k), 



(21) 



(22) 



Pm(k) + P^(k) 



hhr 



(23) 

With the pressure and density field power spectra, one 
can define a bias associated with the large scale pressure, 
relative to density field, 



bu(k)f e = 



> Pn(k) 
Ps(k) 



(24) 



and the dimensionless correlation coefficient between the 
dark matter and baryon distributions 

Pns(k) 



r n (fe) 



(25) 



In Eq. (E 



^P s (k)P n (k) 
the average density weighted temperature 



is 



T, 



dM—^-(M,z)T e (M,z). (26) 
p p am 

Following f3q| , one can define a covariance matrix in 
Fourier space containing the full information on scale de- 
pendence of bias and correlations: 



C(k) = 



( P ss (k) Pus(k) \ _ 



Pus{k) Pnn(k) 



= Pss(k) 



f burn 
burn i>n 



(27) 

The observation measurement of bjj and rn can be con- 
sidered by an inversion of the SZ-SZ, lensing-lensing and 
SZ-lensing power spectra as a function of redshift bins 
in which lensing-lensing or SZ-lensing power spectra are 
constructed. We discuss these possibilities later. 



The calculation of galaxy-pressure power spectrum re- 
quires knowledge on the galaxy distribution within dark 
matter halos. Following, jTHj , we describe the average 
number of galaxies per halo, (N g ) in Eq. (|l2j), such that 



(N g ) = 



M 



0.6 



M > M min 
M < M min 



(31) 



where M„ 



the minimum dark matter halo mass in 



which a galaxy is found, is taken to be 5.3 x 10 n /i _1 
M(7 ) for our fiducial ACDM cosmological model following 
p7fl . The above relation is consistent with semi-analytical 
models, however, we ignore scatter in the observed dis- 
tribution on the mean number of galaxies per halo. 

With the the average number of galaxies per halo, as 
a function of mass, the mean number density of galaxies 
can be written as an integral over the PS mass function 



= JdM{N g )^(M,z). 



(32) 



In practice, the cross-correlation between galaxy distri- 
bution and any other field requires the knowledge on the 
observable galaxy properties such as the magnitude limit, 
relation between luminosity and mass etc. The same 
restriction arising from observing conditions can be in- 
troduced as part of the weight function that takes into 
account the redshift projection of galaxies. 



F. Pressure Bispectrum 
Following p0|], we can write the pressure bispectrum 



R t _ R PPP , R Phh , R hhh 



(33) 



where 

££ pp (fci, k 2 , k 3 ) = ir^l^.l,,. k 3 ) , (34) 
for single halo contributions, 
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k(hMpc" 1 ) k(hMpc" 1 ) 

FIG. 3. Pressure-galaxy cross-correlation power spectra, (a) Comparison of pressure-gas power spectrum to the dark matter 
density-density and galaxy-galaxy power spectra, (b) The variations in pressure power by changing the maximum mass 
considered. The figure follows Fig. hi 



B 



Phh 



(A*, k 2 , k 3 ) = I$f(k u fc 2 )/ 1 ° ff 1 ' (fc 3 )F(fc 3 ) 



Perm. 



(35) 



for double halo contributions, and 



< lh (fcl,/C2,fc 3 ) 



2J(k 1 ,k 2 ,k 3 )lll 1 g '°(k 3 



(36) 

)/ 1 1 „ 1 '°l(fc 2 )P(fc 1 )F(fc 2 )+Pcrm 



schemes to approximate the effect of gas pressure. One 
such scheme that has fractional errors on the 10% level 
for overdensities < 10 is to filter the density perturba- 
tions in the linear regime as P 2 = f^{k/k-p)P 2 s and treat 
the system as collisionless baryonic particles. Their best 
fit is obtained with the filter 



for triple halo contributions. Here the 2 permutations 



are fc 3 <- 
us that 



fci, k,2- Second order perturbation theory tells 



J(ki, k 2 , k 3 ) = 1 



-2/63 



k 2 



k\ 



k 2 



^'1 



k 2 



k\ 



k\ 



2fcifc 2 

2/63 



(37) 



In addition to the pressure bispectrum, one can also 
define cross-correlation bispectra such as the pressure- 
pressure-density or pressure-density-density. These bis- 
pectra are relevant to the calculation of non-Gaussianities 
present in CMB through secondary anisotropies (e.g., 
P0| ) , and is necessary to determine the higher order mo- 
ments associated with cross-correlations between individ- 
ual effects such as SZ and weak lensing. 



G. Baryons in small overdensities 

The power spectrum of baryons that trace the Jeans- 
scale smoothed dark matter density field can be cal- 
culated following (27) (GH), where they studied simple 



fb 



1, 



1 



(1 +4x 2 ) 1/4 



(38) 



and [^J suggests fcp = 34J7m 2 /iMpc _1 as a reasonable 
choice for the thermal history dependent filtering scale. 

For such baryons, we assume that their temperature is 
related to photoionization energy (< 25 eV). The mass 
fraction of baryons present in such small overdensities as 
a function of redshift is obtained through the numerical 
simulations of [Q. 



H. Results &: Discussion 

In Fig. |l](a-b), we show the pressure power power 
spectrum today (z — 0), written such that A 2 (A:) = 
k 3 P(k)/2ir 2 is the power per logarithmic interval in 
wavenumber. In Fig 0(a), we show individual contribu- 
tions from the single and double halo terms and a com- 
parison to the Jeans-scale smoothed dark matter density 
field power spectra, both linear and non-linear, following 
p7[ and using |32j fitting function. Here, we have taken 
the electron temperature to be the virial temperature 
given in Eq. (0). Shown here is also the gas bias bn(k); 
at large scales, bjj(k) ~ 3 as k — > 0, consistent with nu- 
merically measured bias for gas (e.g., jl3]]) and analytical 
estimates in |14|. The pressure power spectrum is such 
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FIG. 4. Pressure bispectra. In (a), the long dashed line is the total density field bispectrum. (b) The variations in pressure 
power by changing the maximum mass considered. An increase in energy does not significantly increase the contribution to 
the bispectrum when halos with high masses are considered. This is due to the strong dependence on rare and most massive 
halos to the bispectrum arising from single halo term. 



that at scales below ~ few h Mpc , the pressure fluctu- 
ations are suppressed relative to the dark matter power 
spectrum; The resulting power spectrum can also be de- 
scribed as a smoothed, but biased, version of the dark 
matter power spectrum. The scale at which smoothing 
enters in to the power spectrum is determined by the 
scale radius of the dark matter and gas profiles. Thus, 
the direct measurement of the pressure power spectrum, 
to some extent, can be used as a probe of halo profiles. 

In the same figure, we also show the measured pressure 
power spectrum in hydrodynamical simulations by fl3| l 
for their ACDM model. For comparative purposes, we 
have appropriately corrected their pressure power spectra 
based on the mean temperature of baryons as tabulated 
p3| since our definition of the pressure power spectrum 
includes the temperature. The resolution of simulations 
limit the accuracy of power spectrum to the range in 
wavenumbers of 0.2 < k < 2.0 h Mpc -1 , and is only 
based on a single realization. In this range, we find that 
our analytical models predict more power than what is 
measured, while agreement is observed at scales of a ~ 
few h Mpc -1 . The extent to which our analytical cal- 
culations agree with simulations is encouraging; this is 
the first time that a detailed analytical model for the 
pressure power spectrum has been compared with a nu- 
merically measured one. Numerical simulations with im- 
proving dynamical range and resolution will eventually 
test the reliability of models such as the one presented 
here as a useful description of the pressure fluctuations 
of the universe. Till then, we consider the present model 
as an appropriate descrption of the large scale pressure 
fluctuations. 



In Fig. [y(b), we show the dependence of pressure power 
as a function of maximum mass used in the calculation 
with maximum mass ranging from 10 16 , 10 15 , 10 14 and 
10 13 Mq from top to bottom. Here, we have shown the 
single halo contribution. Also shown are the total con- 
tribution to pressure power spectrum when there is an 
additional source of energy. Here, we have taken the 
minimum temperature to be ~ 0.75 keV; power spec- 
tra, in general, scale as the square of this energy if the 
real preheating energy is higher or lower than the one 
considered here. There are clear differences between the 
pressure power spectra with and without an additional 
source of energy. With increasing such additional non- 
gravitational energy, note that b in Eq. (||) — > such 
that p g (r) — > p g o- Thus, there is no longer a clear turn 
over in the pressure power spectrum since the effect of 
smoothing resulting from scale radius r s is not present. 
The changes suggest the possibility that physical proper- 
ties associated with large scale structure baryons can be 
probed with pressure power spectrum. In fact, the com- 
bined study of dark matter and pressure power spectra 
may allow a consistent determination of halo properties, 
and to break certain degeneracies associated with dark 
matter halo profile and concentration as noted in [jL5j, 
while at the same time investigating presence of addi- 
tional sources of energy. 

In Fig. U and we study the cross-correlation power 
spectra between pressure and density field and pressure 
and galaxy distribution, respectively. These power spec- 
tra are relevant to the study of correlations present be- 
tween, say, SZ and weak gravitational lensing and SZ 
and galaxies, or a similar tracer of large scale structure, 
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FIG. 5. SZ power spectrum, (a) The halo-halo correlations (dashed line) contributes less than 10% to the total SZ power 
(solid line). The signal halo terms is shown with a dotted line. We also show the SZ power spectra based on our previous 
prediction using a scale-independent bias model of PD non-linear dark matter power spectrum and using the GH Jeans-smoothed 
unbiased power spectrum to describe baryons present in small overdensities with photoionization temperatures. For comparison, 
we also show the power spectrum of lensing primary anisotropies and the upper limits on temperature anisotropies at arcminute 
scales from BIMA and ATCA. (b) The effect of maximum mass on the power spectra, with maximum mass as in Fig. 1. The 
dot-dashed lines are the total SZ power when the minimum temperature is 0.75 keV. We also show the noise variance of the 
Planck SZ map, as has been calculated based on SZ spectral dependance and Planck detector noise. The power spectrum of 
SZ effect will be easily measured with Planck. 



such as radio sources. The presence of additional source 
of energy clearly affects the cross-correlation power spec- 
tra, suggesting the possibility that such effects may be 
investigated using cross-correlations between a tracer of 
pressure fluctuations and a tracer of matter density fluc- 
tuations. 

Since the bispectrum generally scales as the square of 
the power spectrum, it is useful to define 

AcqW^T^W^M), (39) 

which represents equilateral triangle configurations. In 
Fig. |J, we show the pressure bispectrum as produced by 
baryons present in virialized halos. Here, most of the 
contributions at relevant scales come from the single halo 
term. Given the additional dependence on temperature, 
and thus mass, the bispectrum is more strongly sensitive 
to the presence of rare and most massive halos. Thus, an 
increase in energy of such rare halos does not significantly 
change the pressure bispectrum, but such energy changes 
contribute when halos of mass < 10 14 Mq. 

III. SZ EFFECT 

The temperature decrement along the line of sight due 
to SZ effect can be written as the integral of pressure 
along the same line of sight 



y=- =g(x) / dra{r)^—^n e {r)T e {r) (40) 

where ot is the Thomson cross-section, n e is the electron 
number density, r is the comoving distance, and g(x) = 
xcoth(x/2)— 4 with x = hv/ksTcMB is the spectral shape 
of SZ effect. At Rayleigh- Jeans (RJ) part of the CMB, 
g(x) = -2. 

The spectral dependance of SZ effect is unique in 
that it can be separated from most other contributors 
to CMB temperature fluctuations, including the primary 
anisotropy itself. As discussed in Cooray et al. || , the up- 
coming multifrcquency CMB satellite and ballone-borne 
data, among which Planck providing the greatest infor- 
mation on SZ, allow the possibility for detailed studies 
on the SZ effect including its higher order correlations 
such as bispectrum and skewness. Since these observa- 
tions are projected measurements of the pressure power 
spectrum and bispectrum along the line of sight, we now 
provide analytical predictions based for SZ effect based 
on our model for the pressure fluctuations. 

A. SZ Power Spectrum 

The angular power spectrum of the SZ effect is defined 
in terms of the multipole moments yi m of temperature 
fluctuations as 
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FIG. 6. SZ bispectrum. (a) The triple halo, halo-halo-halo, (dot dashed line) and double halo, Poisson-halo-halo, correlations 
contributes less than 10% to the total SZ bispectrum (solid line), (b) The effect of maximum mass on the bispectrum, with 
maximum mass as in Fig. 1. The dot-dashed lines are the total SZ bispectrum when the minimum temperature is 0.75 
keV. Given the strong dependence on mass here, the increase temperature does not significantly change the bispectrum when 
maximum mass is greater than 10 Mq. 
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(41) 



Cf z is numerically equal to the flat-sky power spectrum 
in the flat sky limit. The SZ power spectrum can be 
written as a redshift projection of the pressure power 
spectrum 



C, sz 



dr 



W bZl (rY 



pt 
r nn 



I 



(42) 



where (La is the angular diameter distance. At RJ part 
of the frequency spectrum, the SZ weight function is 



W sz (r) = -2 



kB<JTn e 
a(r) 2 m e c 2 



(43) 



where n e is the mean electron density today. In deriving 
Eq. (42), we have used the Limber approximation El by 
setting k = I /(1a and flat-sky approximation. In previous 
studies (e.g., |ll| and references therein), the SZ power 
spectrum due to massive halos have been calculated fol- 
lowing projected y parameter of individual halos. The 
two approaches are essentially same since the order in 
which the projection is taken does not matter, except 
that our approach allows us to calculate intermediate 3d 
properties of baryons, mainly pressure. 

In Fig. |^(a), we show the SZ power spectrum due to 
baryons present in virialized halos compared with our 
previous prediction for SZ effect using a biased power 
spectrum for pressure fluctuations following |32[ non- 
linear dark matter power spectrum. As shown, most of 
the contributions to SZ power spectrum comes from in- 
dividual massive halos, while the halo-halo correlations 



only contribute at a level of 10% at large angular scales. 
This is contrary to, say, the lensing convergence power 
spectrum discussed in p(J, where most of the power at 
large angular scales is due to the halo-halo correlations. 
The difference can be understood by noting that the SZ 
effect is strongly sensitive to the most massive halos due 
to T oc M 2 / 3 dependence in temperature and to a lesser, 
but somewhat related, extent that its weight function 
increases towards low redshifts. The lensing weight func- 
tion selectively probes the large scale dark matter den- 
sity power spectrum at comoving distances half to that 
of background sources (z ~ 0.2 to 0.5 when sources are 
at a redshift of 1), but has no extra dependence on mass. 
We have also shown current upper limits on the tem- 
perature fluctuations at arcminute scale angular scales 
where potentially the physical properties of baryons can 
be studies. These limits come from Q (BIMA) and E|] 
(ATCA). 

Also shown is the contribution to SZ effect from 
baryons present in overdensities < 10 (curve labeled GH). 
The SZ power spectrum here was calculated by replacing 
the pressure power spectrum in Eq. ( (42| ) with the unbi- 
ased Jeans-scale smoothed dark matter power spectrum 
following |2?J and assuming a mean temperature of 25 eV 
for these baryons. The mass fraction of baryons present 
in such small overdensities were taken from numerical 
simulations of |Q and roughly follows ~ 0.25(1 + z), such 
that at a redshift of 3 and above all of the baryons are 
present in such small overdensities. The power spectrum 
due to such baryons are roughly three orders of mag- 
nitude lower than the power spectrum predicted for SZ 
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FIG. 7. (a) SZ skewness as a function of smoothing scale. The absolute value of skewness, since as 
g(x) — —2 at RJ part of the frequency spectrum, is shown for the virial temperature (solid lines) and minimum temperature 
(dot-dashed lines) models for halos as a function of maximum mass used in the calculation ranging from 10 16 to 10 13 M from 
top to bottom, (b) The probability distribution function when the smoothing scale is 12 arcmin as a function of maximum 
mass used (solid lines; the highest peak curve is when the maximum mass is 10 13 Mq, while with increasing non-Gaussianity as 
demonstrated by the departure from a Gaussian distribution, the maximum mass increases) . We also show the pdf of expected 
Planck SZ map noise (dashed line) for smoothing at same angular scales. The non-Gaussian tail, at the negative y values 
beyond the pdf of Planck noise, due to massive and rare clusters will easily be detected with Planck. 



effect from baryons in virialized halos, but as shown in 
Fig. ||(b), this level is consistent with what is predicted 
for SZ effect when halos with mass greater than 10 13 Mq 
is not present in observed fields. 

As shown in Fig. ||(b), the lack of massive halos leads 
to a strong suppression of power, and halos with masses 
greater than 10 15 Mq are needed to obtain the full power 
spectrum. The lack of massive halos not only lead to a 
change in the power spectrum at large angular scales, the 
lack of mases also affect the contribution at small angular 
scales. Increasing the minimum temperature of electrons 
from the values determined by virial theorem to a min- 
imum energy value of 0.75 keV significantly affects the 
change resulting from the lack of massive halos. In fact, 
with a minimum energy for baryons, the change is smaller 
when halos with masses less than 10 14 Mq is considered. 
At the higher end of masses, the minimum energy do 
not significantly affect the power spectrum; the resulting 
change is less than 30% compared to the power spectrum 
with electron temperature based on the virial theorem. 
The variations suggest several observational possibilities, 
including the determination of minimum electron tem- 
perature, ie. the energy related to preheating if it exists, 
by calculating the power spectrum with massive halos 
substratced in a wide-field SZ map such as the one that 
will be eventually made with Planck. 

For less area surveys, such as planned interferomet- 
ric observations of the wide-field SZ effect (e.g., the few 



square degree survey of Carlstrom et al. |7j]), the sample 
variance due to lack of massive halos in observed fields 
can be problematic in the interpretation of the observed 
signal. The problem arises from the fact that massive ha- 
los which contribute to the SZ power spectrum are rare 
and that observations in small fields will not contain such 
adeqaute masses to provide the fully expected SZ signal. 
The dependance of SZ effect on massive halos is even 
problematic for numerical simulations with limited box 
sizes. As pointed out by 0, the measured power spec- 
trum in their simulation varies significantly based on the 
considred line of sight. 

The dependance of signal on massive halos is also 
present in other observables of large scale structure, such 
as weak gravitational lensing. Compared to weak lensing 
surveys, studied in [^(| and |jl9), the SZ effect depends 
more strongly on rare halos. Most of these halos are 
at low redshifts, thus, surveys which avoid regions with 
known clusters will inherently also include an additional 
bias. As an example, the contribution to 1-a detection 
of temperature anisotropics at arcminute scales by |4|] 
due to SZ effect requires detailed knowledge on the dis- 
tribution of halo masses in the observed fields. For a 
measurement of the SZ power spectrum, with a sample 
variance less than 20%, requires observations of a field ~ 
1000 deg 2 , while the same can be achieved in an area of ~ 
100 deg 2 for lensing. As discussed in poj , however, the 
sample variance due to lack of massive and rare halos, 
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which dominate the SZ power, does not directly imply a 
systematic bias as long as one uses an approach similar to 
the one suggested by carefully accounting for the sample 
variance that may be present from lack of massive halos. 
Such an approach requires a reliable model for the SZ ef- 
fect and detailed numerical simulations will be required 
for such a study. 

These issues can be ignored for the upcoming wide field 
CMB experiments, such as Planck, where the frequency 
coverage will allow the recovery of SZ effect over 65% of 
the sky not confused by galactic emissions, thereby, pro- 
viding an accurate measurement of its power spectrum 
and higher order statistics (see, ||). Such a wide-field SZ 
map is also highly desirable for several reasons including 
the presence of adequate mass distribution of the universe 
such that a fair sample is considered and the possibility 
to use such a wide-field map for various cross-correlations 
purposes, such as against Sloan galax y dis tribu tion o r a 
wide-field weak lensing survey (see, § 1IIC and [II D| ). 



B. SZ Bispectrum & Skewness 



The angular bispectrum of the SZ effect is defined as 

{yhm i yi2m2yi3m 3 } : 

and can be written following Limber approximation as 
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The more familiar flat-sky bispectrum is simply the ex- 
pression in brackets [flij . The basic properties of Wigner- 
3j symbol introduced above can be found in || . 
Similar to the density field bispectrum, we define 



(46) 



we show individual contributions A^ qJ 



involving equilateral triangles in Z-space. The absolute 
value of B sz is considered in above since B sz oc g(x) 3 , 
which is a negative quantity at RJ part of the frequency 
spectrum with g(x) = —2 
In Fig. |(a 

with a maximum mass of 10 le Mq. As shown, the main 
contribution to bispectrum comes from individual halo 
term, while other terms involving correlations between 
halos contribute < 10%. In Fig. [|(b), we show bispec- 
trum as a function of maximum mass used in the cal- 
culation. Here, we have shown the total contribution to 
the bispectrum in solid lines while the total bispectrum 
in the presence of a minimum temperature of 0.75 keV 
is shown wi th a dot-dashed line. The bispectrum, as dis- 
cussed in § II H, is strongly sensitive to the single halo 



term due to additional mass weighing. Almost all of the 
contributions to the SZ bispectrum comes from the single 
halo term. The same dependence on mass massive rare 
halos decreases the effect of a temperature increase when 
maximum mass in > 10 15 Mq. This is in contrast to the 
power spectrum, where differences are still present with 
an increase in temperature from virial to a minimum of 
0.75 keV. 

The measurement of full bispectrum in million to bil- 
lion pixel data of a wide-field SZ map as the one that'll 
be produced with Planck can in general be difficult. In 
fact, there is no algorithm yet to measure the full bis- 
pectrum in such wide-field data in a reasonable time and 
computational requirements. Given such a possibility, it 
is interesting to consider a collapsed measurement of the 
bispectrum; real space statistics such as the third mo- 
ment and skewness allow this possibility. In fact, the 
skewness has now been measured for the COBE data by 
fsof , while the bispectrum measurements have only been 
limited to specific configurations of the bispectrum such 
as equilateral triangles in Z-space The skewness al- 
lows an easily measurable aspect of the bispectrum and 
will be porbably be one of the first measurements of non- 
Gaussianity in a wide-field SZ map. The skewness can be 
calculated using the second, (y 2 (a)), and third, (y 3 (a)), 
moments of the SZ effect: 



(47) 



where the two moments arc 
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In Fig. 0(a), we show the absolute value of skewness, 
I S3 (a) I , as a function of smoothing scale a when the max- 
imum mass included in the calculation ranges from 10 16 
to 10 13 Mq (from top to bottom). The absolute value of 
skewness is considered since £3 oc g(x)~ 1 , which is a neg- 
ative quantity at RJ part of the frequency spectrum with 
g(x) ~ —2. As shown, the SZ skewness is heavily depen- 
dent on the presence of most massive and rare halos. The 
introduction of a minimum energy of 0.75 keV leads to 
decrease in skewness, which results from the fact that the 
power spectrum is more affected than the bispectrum by 
such an increase. 

Given the dependence on rare and most massive halos, 
the SZ effect is considerably non-Gaussian. As discussed 



12 




Gaussianity, one can write the pdf of SZ effect to second 
order as 



FIG. 9. The cross-correlation coefficient for SZ-weak lens- 
ing (thick lines) and SZ-galaxy (thin lines) with (solid line) 
and without additional energy (dashed line). In general, when 
the maximum mass is greater than 10 14 Mq, SZ and galaxy 
correlate less than SZ-weak lensing at small angular scales, 
while the opposite is present at large angular scales. Introduc- 
tion of additional energy increases the correlation by not more 
than few percent; such an increase is insignificant to be de- 
termined observationally using the SZ-lensing and SZ-galaxy 
cross-correlations, however, SZ power spectrum and skew- 
ness allow a determination as these values are significantly 
affected. 



in Pq| , the non-Gaussianity can be used as a useful tool 
for the identification and separation of most rare and 
massive halos from the wide-field CMB data. An opti- 
mised algorithm that utilised both the non-Gaussianity 
of SZ effect and its frequency dependance will be use- 
ful for constructing a catalog of SZ clusters in upcoming 
wide-field data. With massive clusters separated out, the 
remaining contribution to SZ effect will be from halos of 
mass < 10 14 Mq, such as galaxy groups. In we con- 
sidered contribution from such small halos as the one due 
to large scale structure. The extent to which such small 
halos contribute to the SZ power spectrum and higher 
order statistics clearly depend on the role of additional 
energy in baryons. Therefore, any measurement of the 
power spectrum with known massive clusters removed, 
can be in return used as a probe of physical properties 
related to baryons, mainly the extent to which prehating 
affects the electron temperature of low mass halos. 

Instead of individual non-Gaussian statistics as skew- 
ness, one can construct the probability distribution func- 
tion (pdf) using a SZ map smoothed on some scale a. 
The use of pdf as a probe of cosmology was first suggested 
by fl5| for weak gravitational lensing convergence. The 
same technique can be easily extended to SZ. Using the 
Edgeworth expansion to capture small deviations from 
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where H 3 (x) = x 3 — 3x is the third order Hermite poly- 
nomial (see, Q for details). 

In Fig. 0(b), we show the pdf of SZ effect at 12' as a 
function of maximum mass used in the calculation. As 
shown, the greatest departure from non-Gaussianity oc- 
cur when the maximum mass of halos are greater than 
10 14 Mq. Given that we have only constructed the pdf 
using terms only out to skewness, the presented pdfs 
should only be consider as approximate; With increas- 
ing non-Gaussianity behavior, the approximated pdfs are 
likely to depart from true distributions especially in the 
tails. Observationally, the pdf can be constructed easily 
by considering a histogram of the pixel temperature val- 
ues of the SZ map. Though such a construction sounds 
straightforward, there is likely to be complications on the 
interpretation of such a pdf in the presence of instrumen- 
tal noise and other foregrounds. Techniques which do no 
directly make a wide-field map, especially interferomet- 
ric observations, will again require special techniques to 
construct the pdf. Therefore, the extent to which the full 
pdf, or such an histogram, can be used a probe of cosmol- 
ogy and the accuracy to which pdfs can be constructed 
from upcoming wide-field CMB anisotropy data, such as 
Planck and planned interferometric surveys (Carlstrom, 
private communication), need to be investigated in de- 
tail. We leave these issues for further study. 



C. SZ-Weak Lensing Cross-correlation 

Similar to the SZ power spectrum, the angular power 
spectrum of weak lensing convergence can be defined in 
terms of the multipole moments ni m as 



(51) 



and can be written in terms of the dark matter power 
spectrum by 



I 

d A ' 



(52) 



When all background sources are at a distance of r s , the 
lensing weight function becomes 



era 



d,A{r s 



(53) 



The detail properties of lensing statistics, under the dark 
matter halo approach, is discussed in and po| . 
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FIG. 8. SZ-weak lensing cross-correlation (a) and SZ-galaxy cross-correlation (b) power spectra as a function of maximum 
mass, with mass cut off following Fig. 1. The additional minimum energy, shown with a dot-dashed line, leads to an increase 
in the correlated power. 



The cross-correlation between the SZ effect and weak 
gravitational lensing can be similarly defined in terms of 
the individual multipole moments as 



( K lm,yi'm') — Cf Z K Sii'5 mm i 



(54) 



This is now related to the dark matter-pressure power 
spectrum by 



1 1 dr A m W 



(55) 



Finally the cross-correlation coefficient between SZ and 
weak lensing is 



C(SZ,ac) ; 



c 



sz- 



cf z c? 



(56) 



D. SZ-Galaxy Cross-correlation 

Similar to the SZ-weak lensing cross-correlation, one 
can study the cross-correlation between the galaxy dis- 
tribution, which traces the large scale structure, and the 
SZ effect. The power spectrum of galaxy distribution can 
be defined terms of the multipole moments gi m as 



{gimdl'm') = Cf5 w 6 n 



(57) 



and can be written as a projection of the dark matter 
power spectrum 
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d A 



(58) 



Here, Cf should be understood as the 2d Fourier trans- 
form of the galaxy correlation function, generally referred 
to as w{9) in the literature. The weight function for 
galaxy projection involves the redshift distribution 

dN 

W K {r) = V (59) 
dr 

normalized such that J dr(dN/dr) = 1. 

As before, the cross-correlation between the SZ effect 
and galaxy distribution can be similarly defined in terms 
of the individual multipole moments as 



(gimVl' 
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(60) 



This is now related to the dark matter-pressure power 
spectrum by 
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Finally the cross-correlation coefficient between SZ and 
weak lensing is 

C(SZ, 5 ) ; = ^^= (62) 



Cf*C? 

In Fig. ^, we show the SZ-weak lensing and SZ-galaxy 
cross-correlation power spectra as a function of maximum 
mass used in the calculation, while the correlation coeffi- 
cients are shown in Fig. ^| In order to describe the galaxy 
distribution, we have considered a survey at low redshifts, 
similar to the Sloan Digital Sky Survey (SDSS)0. Such a 
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low redshift tracer is desirable since contributions to SZ 
effect primarily comes from large scale structures at red- 
shifts < 1. Here, we have assumed the redshift distribu- 
tion of Sloan galaxies follow dN/dz oc z 2 exp[(— z/z c ) 3 / 2 ] 
with a mean redshift z m , of 0.2 (z c ~ 1.412z m ). For SZ- 
weak lensing, the cross-correlation is such that SZ and 
lcnsing traces each other out to angular scales of ~ 1000 
when most massive and rarest halos are involved with a 
decrease in cross-correlation between the two at small 
angular scales. For SZ-galaxy cross-correlation, there 
is additional correlation at large angular scales, when 
compared to lensing, while the correlation is suppressed 
at small angular scales. The decrease at small angular 
scales is due to the fact that small halos that contribute 
to SZ and lensing do not necessarily contribute to the 
galaxy power spectrum. With an increase in additional 
energy for halos, the cross-correlation increases by few 
percent, however, this small increase unlikely to be de- 
termined accurately through observations. The correla- 
tion is sensitive to the redshift distribution of galaxies, 
which depends mostly on the selection criteria imposed 
by observations. The selection function of weak lensing 
can be considered well understood, however, a straight 
forward interpretation of any observed SZ-galaxy cross- 
correlation will require a clear understanding of observ- 
able related to galaxy distribution. 

The SZ-SZ, lensing-lensing and SZ-lensing power spec- 
tra dependent different on the bias and correlation pa- 
rameters. Since the bias and correlation are scale and 
redshift dependent, the measurement of these power 
spectra, which are projected along the redshift distribu- 
tions, do not allow a direct probe of these quantities. 
A useful approach would be to consider the inversion of 
these power spectra in redshift bins by considering the 
measured lensing-lensing and SZ-lensing power spectra 
as a function of redshift. Note that the SZ-SZ power 
spectrum cannot be easily separated in redshift space as 
we do not have the ability to separate individual redshift 
contributions, unlike say in lensing, where one can use the 
redshifts of background sources to construct convergence 
as a function of redshift. With adequate signal-to-noise 
from wide-field surveys, it is likely that such an approach 
will allow studies to be carried out on the extent to which 
temperature weight baryons trace the dark matter and 
their correlation properties. 

Similarly, as studied in |l9]], the cross-correlation of SZ 
against galaxy data, as a function of redshift, is expected 
to provide information on the properties of clustering of 
galaxies with respect to the temperature weighted baryon 
field represented by SZ effect. Such a cross-correlation 
will help understand the extend to which hot /warm gas 
is present in the outskirts of individual galaxies. With 
individual sets of SZ, weak lensing and galaxy maps, it 
is likely that a tremendous amount of information on 
physical properties associated with dark matter, baryon 
and galaxy distribution will be obtained through both 
a comparison of individual power spectra and higher or- 
der moments and cross-correlations and higher order mo- 



ments associated with such cross-correlations. We hope 
to study some of these possibilities in detail in future 
studies. 



IV. SUMMARY & CONCLUSIONS 

Using an extension of the dark matter halo approach, 
we have presented an efficient method to calculate the 
large scale structure pressure power spectrum and its 
high order moments, such as bispectrum. We have di- 
vided the contribution to large scale pressure power spec- 
trum based on the overdensities in which contributing 
baryons are present with (1) baryons present in virial- 
ized halos with overdensities greater than ~ 200 and in 
hydrostatic equilibrium with the density field of such ha- 
los, (2) photoionized baryons in overdensities less than ~ 
10 and which trace the Jeans-scale smoothed dark matter 
density field, and (3) baryons present in the mid overden- 
sity regime which are likely to be undergoing collapse and 
shock heating. 

Our approach allows us to calculate not only 2d statis- 
tics such as the projected pressure power spectrum, or the 
SZ effect, which will be observed, but also the 3d statis- 
tics that can be directly compared to predictions based 
on numerical simulations. We have performed such a 
comparison to recently published numerical simulations 
by |]13| and found good agreement between our analytical 
calculations and their simulations. The current simula- 
tions are limited to a handful of realizations and limited 
dynamical range and resoltion. With improving reso- 
lution and accuracy, analytical models such as the one 
presented here will be tested in detail against numerical 
calculations. Analytical calculations, aided by numerical 
simulations, will eventually allow detailed studies of large 
scale baryon distribution using observations such as the 
wide-field SZ effect. 

The projected pressure power spectrum along the line 
of sight, provides a direct calculation of the large scale 
structure SZ effect and its higher order correlations. In 
the absence of massive and rare halos, we have suggested 
that baryons present in small overdensities provide a 
lower limit to any contribution to SZ effect. The extent 
to which baryons present in overdensities between 10 and 
200 contribute to the correlations in large scale pressure 
and, from it, the SZ effect, requires additional studies, 
preferably with numerical simulations. Presently, our 
Understanding the role of preheating and its effect of 
baryons will also be another challenge as the SZ observa- 
tions will clearly depend on such additional energy con- 
tributions to large scale baryon distribution. We have 
suggested the possibility of using SZ power spectrum and 
higher order correlations, such as the SZ skewness, as a 
probe of preaheating. Such a study will require a wide- 
field SZ map and this task will be completed with Planck 
observations. The unique frequency dependance of the 
SZ effect, together with its non-Gaussian behavior, will 
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allow the construction of a reliable SZ cluster catalog 
with will aid in cosmologicalo studies of structure forma- 
tion. 

Our approach allows one to study possible systematic 
effects that may be present in upcoming SZ observations 
of small area fields due to the presence or absence of rare 
massive halos in such fields that will be observed. We 
have shown that the SZ effect as well as its non-Gaussian 
properties are mainly due to the most massive and rarest 
virialized halos in the universe. The lack of massive halos 
in observed SZ fields can introduce a systematic bias in 
the power spectrum, but the sample variance introduced 
by the lack of such masses, can be easily corrected based 
on the prior knowledge on mass distribution of observed 
fields. Due to additional mass dependence through tem- 
perature, the effect of mass is such that the SZ effect is 
more dependent on the rare halos than weak gravitational 
lensing convergence. The same SZ halos also contribute 
to lensing convergence and the cross-correlation between 
SZ and lensing can be used a probe of clustering prop- 
erties between density and temperature weighted baryon 
fields. Given the great potential to study baryon distri- 
bution using SZ, various issues suggested here involving 
such correlations merit further study. 
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